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We study the spectrum of massive excitations of the three-dimensional <f> 4 and Ising models, in the broken- 
symmetry phase. Using a variational method, we show that the spectrum contains all the + states that one 
expects from duality with the glueball spectrum of the 2Z2 gauge model. From the point of view of continuum </> 4 
theory, we show that at least one of the states we find has a non-perturbative origin. 



1. DUALITY AND UNIVERSALITY 

The three-dimensional Ising model is related 
to two other important three-dimensional theo- 
ries: by duality, to the 2Z2 gauge model, and by 
universality, to <j) 4 theory. 

Duality is an exact equality of partition func- 
tions: it means that the Ising model and the 
7Z>2 gauge model are different descriptions of the 
same physics. In particular, the broken symme- 
try phase of the Ising model is equivalent to the 
confined phase of the gauge model. 

Universality tells us that the Ising model and 
(j> 4 theory behave in the same way when approach- 
ing the critical point. Universal quantities are 
the same in a critical region around the transi- 
tion point. Indeed, the universal quantities of 
the Ising universality class have been predicted to 
great accuracy using perturbative 4> 4 theory (see 
e.g. Ref. §). 

Let us apply the tools of universality and dual- 
ity to the problem of determining the spectrum of 
massive excitations of the 3D Ising model in the 
broken symmetry phase. The glueball spectrum 
in the Zi gauge model has been thoroughly stud- 
ied numerically: it turns out to be a rich spectrum 
with many excitations in various angular momen- 



tum channels. On the other hand, we certainly 
do not expect to find an interesting spectrum in 
perturbative <j) theory, which describes just one 
particle. 

Therefore it seems that duality and universal- 
ity lead to contradictory expectations about the 
spectrum of the Ising model and <p A theory in the 
broken phase. This work clarifies these issues by 
a numerical evaluation of the spectrum of both 
models, performed with a new variational proce- 
dure. 

2. MONTE CARLO DETERMINATION 
OF THE SPECTRUM 

The spectrum of a model is extracted from 
Monte Carlo simulations by studying the long dis- 
tance exponential decay of correlation functions. 
It is convenient to study time slice observables: 
for example if <f> is the order parameter one de- 
fines 
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where L is the lattice size in the x and y direc- 
tions. One then studies connected correlators of 
the S operator, the advantage being that they 
behave as a pure exponential in the long distance 
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limit: 



(S(0)S(t)) c 
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For a theory with a non-trivial spectrum one ex- 
pects to see subleading exponentials as well: 
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Therefore a procedure to determine the spectrum 
from Monte Carlo simulation is to measure time 
slice correlations and fit them to Eq. (j|). 

However there is a more effective method, in- 
spired by what is commonly done in lattice gauge 
theory to determine the glueball spectrum |2j]. 
One introduces a basis of suitably defined (time 
slice) operators Oi (t) and then computes the ma- 
trix of cross-correlators 



G y (t) = <4(0)O, (i)) c 
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Cij(t) is then diagonalized to read off the spec- 
trum. The crucial point is of course the choice of 
the operator basis, which must be carefully fine- 
tuned to obtain an efficient determination of the 
spectrum. Our choice is described in detail in 
Ref. ||. Here we just mention that we included 
the standard magnetization Eq. (Q) in the basis 
{Oi} and that the other operators are defined on 
different length scales. 

We simulated both the Ising model and the lat- 
tice regularized 4> 4 theory in the broken symme- 
try phase, at various temperatures well inside the 
scaling region where universality is expected to 
hold. We considered the + channel only, ft 
turns out that three states can be identified in 
this channel. Scaling is perfectly statisfied since 
the ratios between the masses of the three states 
do not change with the temperature within the 
scaling region. Universality is satisfied as well 
since we obtain compatible results for the ratios 
from the Ising and the 4 simulations. Therefore 
we can quote a single result for each mass ratio: 
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= 2.45(10) 
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Also duality is satisfied since in the TL-i gauge 
model one obtains mass ratios of 1.88(2) and 



2.59(4) in the + channel of the glueball spec- 
trum g. 

Note that the first excited state lies below the 
pair production threshold: this means that, in 
terms of continuum 4 theory, it cannot be of 
perturbative origin. On the other hand the state 
at 2.45 times the fundamental mass could well be 
a signature of the cut in the Fourier transform of 
the propagator induced by self interaction effects, 
and therefore a perturbative effect. 

Let us analyze in more detail what we expect 
from perturbative 4 theory in this respect: a 
simple one-loop computation (see Ref. ||) shows 
that the cut in the momentum space propagator 
implies for the time slice correlators the behavior 

e -2mi|t| 



{S(0)S(t)) c 



-uii \t\ 



c 2 - 



f 



(7) 



However the second term can be shown to be nu- 
merically indistinguishable from an exponential 
decay with mass ~ 2.4 mi. Therefore we con- 
clude that while the state m 3 could be explained 
as a perturbative self- interaction effect, the state 
m,2 is certainly of non-perturbative origin. 

3. BACK TO THE SPIN-SPIN CORRE- 
LATOR 

In this section we apply the knowledge we have 
gained of the spectrum of the theory to the anal- 
ysis of the (time slice) spin-spin correlator 

G(t) = (S(0)S(t)) c (8) 

It is useful to define the effective correlation 
length 

l/£ e// (t)=logG(i)-logG(i + l) (9) 

so that for t — > oo £ e // - * l/"i, Tti being the 
fundamental mass. Clearly if G(t) had a purely 
exponential behavior then £, e ff{t) would be con- 
stant: the preasymptotic behavior of £, e ff(t) de- 
pends on higher mass states and/or interaction 
effects. These effects are shown in Fig.l, where 
m £ e / / is shown for various values of the temper- 
ature in the Ising model. The data from the var- 
ious temperatures are perfectly compatible with 
each other, signaling that the preasymptotic be- 
havior of £ e ff(t) is a physical, scaling effect and 
not a lattice artifact. 
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In Fig. 2 we use our knowledge of the spectrum 
to describe the behavior of tj, e ff(t): the Monte 
Carlo are compared to the perturbative predic- 
tion 

&//(*) =ci [e- mit + f cut (t)] (10) 

(dotted line), and to the curve 

&//(t) = Cl [e- mit + / c „ t (i)] + C2 e- m2t (11) 

(solid line) where the constatnts mi, ni2, ci, c 2 
are taken from our variational evaluation of the 
spectrum, and f cu t is taken from one-loop pertur- 
bative calculations in the continuum theory (see 
Ref. ||). The good agreement between this curve 
and the data suggests that in fact the third mass 
m3 is not a new state but a perturbative interac- 
tion effect. 



4. CONCLUSIONS 

The main result of our analysis is that 3D </> 4 
theory has a rich spectrum of massive excitations 
that signals the existence of non-perturbative 



physics. This spectrum matches accurately the 
corresponding spectrum of the 31? Ising model, 
to which 4 is related by universality, and the 
glueball spectrum of the 3D Z5 2 gauge model, re- 
lated by duality to the Ising model. 

We are currently investigating higher spin ex- 
cited states, corresponding to higher spin glue- 
balls, and the effect of non-perturbative physics 
on the field theoretic prediction of universal quan- 
tities. Another interesting development would be 
to investigate the same issues in other 3D univer- 
sality classes, in particular in TV-component cj) A 
theory. 
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